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Abstract — In textbooks on electromagnetism, the plane wave 
solutions of the Maxwell equations in free space are in most cases 
limited to homogeneous waves where the field is uniform in any 
plane perpendicular to the propagation of the phase. However, 
non homogeneous plane waves where the field is evanescent in a 
direction perpendicular to the propagation also satisfy the 
Maxwell equations. This paper revisits such waves that 
significantly differ from the homogeneous waves. In particular the 
electromagnetic field is not transverse, it has a component in the 
direction of propagation of the phase. The equivalents of the Snell 
law and of the Fresnel coefficients are derived for the reflection 
and transmission of evanescent plane waves at interfaces, and it is 
Shown that such waves can be generated in the computational 
domain of numerical methods that solve the Maxwell equations. 


Index Terms—Electromagnetic theory, Electromagnetic fields, 
Electromagnetic propagation, Inhomogeneous waves, Evanescent 
waves, Plane waves, Maxwell equations. 


I. INTRODUCTION 


The homogeneous plane wave solutions of the Maxwell 
equations in a vacuum or in other infinite media, also called 
traveling waves, play a fundamental role in electromagnetism. 
The simplicity of the solutions permits theoretical derivations 
to be performed, concerning the propagation, the reflection, and 
the refraction of the electromagnetic energy. The plane wave 
concept is thus essential both for the understanding of most 
electromagnetic phenomena and in engineering sciences 
relying on electromagnetism. 

Beside the traveling waves, non homogeneous plane waves 
are also solutions of the Maxwell equations in free space. Such 
waves widely differ from traveling waves. The electromagnetic 
field is not transverse, it possesses a component parallel to the 
direction of propagation of the phase, and it 1s not uniform in 
the planes perpendicular to the direction of propagation, where 
it varies exponentially. 

The non homogeneous plane waves, also called evanescent 
plane waves, are rarely addressed in textbooks on 
electromagnetism, even at research level. Especially, the 
existence of these waves is not mentioned in [1, 2, 3, 4]. The 
only textbook we have found on this question is the one of G.S. 
Smith [5], where a chapter is devoted to non homogeneous 
waves. Concerning the literature, we have not found any paper 


describing with details such waves. We published a paper in 
1999 [6] where evanescent plane wave solutions are derived in 
a vacuum and in the Perfectly Matched Layer (PML) medium. 
However, [6] is limited to the special two-dimensional (2D) 
case that allows the reflection observed from the PML 
absorbing boundary condition to be interpreted. 

In the present paper, we revisit the evanescent plane waves 
in the general three-dimensional case and we derive the 
corresponding Snell Law and Fresnel reflection and 
transmission coefficients at interfaces between dielectric media. 
We also show that evanescent waves can be generated in the 
computational domain of numerical methods, which is of 
interest to assess theoretical derivations and to perform 
experiments involving evanescent fields. 

More specifically, in section II the evanescent plane wave 
solutions of the Maxwell equations are derived in a manner 
which is different, and we think simpler, than that in [5]. The 
waves are described with details, in particular the two possible 
modes where either the electric field or the magnetic field has a 
component parallel to the propagation of the phase. Explicit 
expressions of the wavenumbers and of the field components 
are provided in Cartesian coordinates, which allows such 
problems as the transmission and reflection of evanescent 
waves at interfaces to be easily addressed. 

In section III the Poynting vector and the relationships that 
connect the EF, H, and k vectors of the waves are derived. 
Section IV is devoted to numerical experiments with the finite- 
difference time-domain (FDTD) method [7]. It is demonstrated 
that evanescent plane waves can be enforced in the FDTD grid 
where they propagate in accordance with theory. In section V, 
it is shown that the evanescent fields in some problems, in 
particular in waveguides and on corrugated surfaces, can be 
interpreted as a set of evanescent plane waves. 

Finally, sections VI and VII address the question of the 
reflection and transmission of evanescent plane waves at the 
interface between dielectric media. We have not found any 
textbook or paper on this question, except for the special case 
where an evanescent wave is present behind the interface when 
the incident traveling wave is reflected in totality. Here, the 
Snell law and the Fresnel coefficients are derived in the general 
case where the incident plane wave is evanescent. The reflected 
and transmitted fields are discussed with details, in particular in 
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some special cases. The theoretical Snell law and Fresnel 
coefficients are verified by means of numerical experiments 
with the FDTD method in the case where the evanescence is 
lying in the plane of incidence. 


Il. EVANESCENT PLANE WAVE SOLUTIONS OF THE MAXWELL 
EQUATIONS IN FREE SPACE 


In this section the waveform and the structure of the field of 
the plane wave solutions of the Maxwell equations are derived 
and discussed. 


A-The wavenumbers and the waveform 


Let us consider harmonic plane waves with the 
electromagnetic field assumed in the form 
E jot —jkr 

y(r t) =e (2-1) 


where y is either the E vector, the H vector, or any component 


of E or H, and where is the angular frequency, k is the wave 


vector, and F is the position vector. In coordinates (x, y, z) the 
waveform reads 

jot—jk,x-jk,, y—jk, z 

W(x, V,Z,t) =Wye i (2-2) 

where kx, ky, kz- may be real or complex components of the wave 

vector. By enforcing (2-2) into the Maxwell equations, the 


equation of dispersion is obtained as: 
2 


Kak pke (2-3) 
C 

where c is the speed of light. For waveform (2-1)-(2-2) to be a 
solution of the Maxwell equations, wavenumbers kx, ky, kz must 
satisfy (2-3). 

Equations similar to (2-2) and (2-3) also hold in any other 
coordinates, in particular in coordinates (X, U, V) in Fig. 2-1 
with X axis in direction (0, ọ) 


2 
k? +k +k? =Z (2-5) 
X U V c? 


where kx, ku, ky are connected to kx, ky, kz- by the change of 
coordinates from (x, y, z) to (X, U, V). 

Let us search solutions (kx, ku, kv) that satisfy (2-5). This will 
provide us with plane wave solutions (2.4) of the Maxwell 
Equations. A trivial solution to (2-5) 1s 


k,=t—; ky=0; k,=0 (2-6) 
Choosing the positive sign of kx yields the waveform 
jøt-jÊX 
W(X,U,V,th=we i (2-7) 


This is a traveling plane wave whose phase propagates in 
direction X (-X if the negative sign of ky were chosen) and 
magnitude is uniform in planes perpendicular to _X. Using (2-7) 
into the Maxwell equations expressed in (X, U, V) coordinates, 
it is easily shown that the components of F and H fields are 
perpendicular to the direction of propagation X and are 
perpendicular to each other. 





Fig. 2-1. The system of coordinates (X, U, V) with X in direction (9, ọ). 


Let us consider again the equation of dispersion (2-5). 
Solution (2-6)-(2-7) is not the sole solution. It is trivial to verify 
that the following solution also satisfies (2-5) 


ky = +2 cosh y (2-8a) 
c 
O. 
k, =+j—sinh y cosy (2-8b) 
c 
k, =+ j— sinh y sny (2-8c) 
C 


where y and 1 are two arbitrary real numbers, and where the 
arbitrary signs of ky and ky should be identical and are 
independent of the sign of ky. Choosing the positive sign for ky 
and the negative sign for ky and ky, and inserting these 
wavenumbers into (2-4), the following waveform is obtained 


-j2 coshy X se sinhy(cos7U+sinnV ) 
c e 


W(X, U, V) = Wo e (2-9) 





X 


Fig. 2-2. The systems of coordinates (X, U, V) and (X, Y, Z). 


Defining now direction Y lying in the (U, V) plane and forming 
the angle n with U axis (Fig. 2-2), we have cos7 U+ sinn V = 
Y, so that (2-9) becomes 


Ege cosh y X 


W(X,Y,Z)=yw,e * e (2-10) 


The wavenumbers (2-8) can then be rewritten in the (X, Y, Z) 
Cartesian coordinates as 


k, =“cosh y: k, =—j “sinh 7; k,=0 (2-11) 
c c 


a sinhy Y 
c 


The plane wave solution of the Maxwell equations in free 
space (2-10) is thus a wave whose phase propagates 1n direction 
X with celerity c/coshy, and whose magnitude is constant in 
direction Z and evanescent in direction Y, where X, Y, Z are 
orthogonal. Obviously, by means of other choices of the 
arbitrary signs in (2-8), the signs in (2-10)-(2-11) could be 
changed and then the wave could propagate in direction —X¥ or 
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be evanescent in direction —Y. We can notice that changing the 
sign of y is equivalent to changing the signs of (2-8b) and (2- 
8c). In this paper we assume that the signs are those in (2-10)- 
(2-11), with y defined as positive or negative which permits the 
direction of evanescence to be either + Y or —Y for a direction of 
propagation assumed in +X direction. 

For a specified direction of propagation X, the angle n in (2- 
8) is arbitrary so that the direction of evanescence Y can be any 
direction perpendicular to X. The evanescence is characterized 
by the parameter y. For y = 0 the magnitude is uniform, the 
wave is a pure traveling wave, 1.e. (2-11) reduces to (2-6). For 
X # 0, the larger the modulus of y the stronger the evanescence. 


t/2—0 
z A | 





Fig. 2-3. The directions X of propagation of the phase, Y of evanescence, and 
Z. Starting from axes (x, y, z) a rotation by ~ around z gives axes (x’, y’, z), from 
which a rotation by 7/2- 8 around y’ gives axes (X, y’, z’), and finally a rotation 
by n around_X gives orthogonal axes (X, Y, Z). 


Waveform (2-10) and wavenumbers (2-11) express the 
evanescent wave in its own system of orthogonal coordinates, 
where X and Y are the directions of propagation and 
evanescence, respectively. Obviously, it 1s possible to express 
a given physical evanescent wave in any other system. For 
instance in the system (x, y, z) represented in Fig. 2-3, where 
propagation X is in direction (8,@), and evanescence in 
direction Y defined with the angle n. By applying the change of 
coordinates from (X, Y, Z) to (x, y, z) to the wave vector, that is 
to wavenumbers (2-11), wavenumbers kx, ky, k, are obtained as: 
k = 2 Teosh y cos osin 0 

C 


X 


+ j sinh y (cos o cos 0 sin y + sin p cos n)| (2-12a) 


w er 
k, = — [cosh z sin psin 6 
+ j sinh y (sin g cos 0 sin y — cos Q cos n )| (2-12b) 
k, = “Icosh y cos 0 — sinh ysin Osin n] (2-12c) 
c 


It can be verified that the set (2-12) satisfies the equation of 
dispersion (2-3). And for y = 0 it reduces to the wavenumbers 
(@ cos@ sinO/c, w sing sin8/c, @ cos@/c) of a pure traveling 
wave propagating in direction (9, ©). 


B-The field components of evanescent plane waves 


With uniform plane waves, £ and H are perpendicular to each 
other and perpendicular to the direction of propagation. This is 
no longer true with evanescent waves. To show this, let us 
consider the Maxwell equations expressed in the (X, Y, Z) 
coordinates. And let us enforce the six components of the field 
in the form (2-10) multiplied with exp(jœt), with component 
magnitudes Exo, Evo, Ezo, Hxo, Hyo, Hzo. The following S1X 
equations are obtained 


ECE nm Jsi ZH (2-13a) 
E,cE,, =coshy H,, (2-13b) 
é,cE,,=-coshy H,,—jsinhy H,, (2-13c) 
lo CH y =—j sinh y Ezo (2-13d) 
La CHp =—cosh y E,, (2-13e) 
Lo CH7 = cosh y Eyo + j sinh y Eyo (2-13f) 


We can note that (2-13) is composed with two independent 
subsets, (2-13a)-(2-13b)-(2-13f) and (2-13c)-(2-13d)-(2-13e), 
respectively. This means that two independent solutions can 
propagate in direction X and be evanescent in direction Y. 
Subset (2-13a)-(2-13b)-(2-13f) yields 


E = jZ sinb y Ho (2-14a) 
Ezo =Z cosh y H, (2-14b) 
Hy =H, (2-14c) 
and subset (2-13c)-(2-13d)-(2-13e) 
EE. (2-15a) 
Po Sing (2-15b) 
Zo 
,, = > cosh y E, (2-15c) 


0 

where Ho and Æo are arbitrary constants and Zo is the impedance 
of free space. We denote the two solutions as TM mode and TE 
mode, with respect to the direction of propagation, respectively. 

The TM mode (2-14) has Ey and Az components 
perpendicular to the direction of propagation, as uniform plane 
waves. But the ratio of the components is different, Z) coshy in 
place of Zo. There is an additional Ey component in the direction 
of propagation. Due to the presence of j = exp(jm/2) in (2-14a), 
Ey is in advance by 7/2 with respect to Ey. It vanishes for y =0 
and equals Ey in magnitude when y is large. The physical 
components corresponding to (2-14) are the real parts of the 
product of (2-14) with exp(-j@ coshy X/c), exp(-w sinhy Y/c), 
and exp(jaf). This yields 


E,(X,Y,0 =-Z, H,sinhy A, sin(at—¢, ) (2-16a) 
Ey (X,Y,t)=Z, H, coshy A, cos(at—¢, ) (2-16b) 
H (X,Y, t) =H, A, cos(at—¢, ) (2-16c) 
where 
A, =e marie 5g = © cosh y X (2-16d) 
C 
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In the TE mode (2-15), the E field is perpendicular to the 
propagation and the evanescence. There is a Hy component in 
the direction of evanescence, and a Hy component in the 
direction of propagation, in advance by 7/2 with respect to Hy. 
This longitudinal Hy vanishes when y = 0. The physical 
components corresponding to (2-15) read 


E,(X,Y,t) =E, A, cos (at -—¢,) (2-17a) 
H,(X,Y,0) =E,/Z, sinh y A, sin (at — ¢, ) (2-17b) 
H, (X,Y, t)=— E, /Z, cosh y A, cos (at — gy ) (2-17c) 


The TM Mode 
i The TE Mode 
Z^ 
i 1 A 
A y Z'ı 
Hz i 
Ey Ez I 
a Y 
--- > if 
Ex X San Za a > 
Hx X 
Hy 


Fig. 2-4. The two evanescent modes with propagation of the phase in X 
direction and evanescence in Y direction. The physical components (2-16)-(2- 
17) are represented for -7/4 < œt-px < 0, i.e. cos(œt-py) > 0 and sin(@t-ox) < 0. 


The field components of the two evanescent modes are 
represented in Fig. 2-4, with either £ or H in the direction of 
propagation. The two modes are dual which results from the 
duality of the Maxwell equations. When y vanishes they 
degenerate into the single traveling wave mode. For each mode, 
the direction of evanescence is any arbitrary direction Y 
perpendicular to the direction of propagation X. The general 
solution of the Maxwell Equations is any linear combination of 
the two modes. 

The field components (2-14)-(2-15) can be expressed in any 
other coordinates, especially in the (x, y, z) coordinates of Fig. 
2-3 where (X, Y, Z) are defined with 0, , n. For the TM mode 
(2-14), the magnitudes of E and H components read 


Ej” |- (cos pcos Osin7 + sin pcos 7) cosh y 
+ jcosp sin@ sinh y|Z,H, (2-18a) 


Ey = (cos pcos —sing cos @sin7)cosh y 


+ jsing sin@sinh y|Z,H, (2-18b) 
E_, =|sin@sinn cosh y+ jcos@sinh y|Z,H,  (2-18c) 
H „ = [sin gsinņ — cos gcos 6 cos n| H, (2-18d) 
H „ = -cos osin y + sin pcos 6 cos 7) H, (2-18e) 
H =sinĝcosy H, (2-18f) 


where Ho is the magnetic component perpendicular to the 
direction of propagation and the direction of evanescence. For 
the TE mode, (2-15) yields 


E „ = [sin osin y — cos g cos 8 cos 7) |E, (2-19a) 
Evy = —|cos osin y + sin ocos 0 cos nlE, (2-19b) 
E =sinĝ cosy Ep (2-19c) 
A= l(sin pcos n + cos ocos Osin n) cosh y 
— jcos osin @sinh y |E, /Z, (2-19d) 
H o = —|(cos pcos y — sin gcos 8sin77) cosh y 
+ jsingsin@sinh y|E,/Z,  (2-19e) 
H_,,, =—|sin@sinncosh y + jcos@sinh y|E,/Z, 2-199 


where £ 1s the electric field perpendicular to the direction of 
propagation and the direction of evanescence. 

The expressions of the field components (2-18) and (2-19) 
are used in section VII for deriving the reflection and 
transmission coefficients of the evanescent plane waves at 
interfaces between dielectric media (Fresnel coefficients). 

It can be noticed that modes (2-14) and (2-15) are simpler in 
appearance than the modes in the textbook [5]. This is because 
they are expressed in the natural coordinates (X, Y, Z) of the 
waves. It can be shown that (2-14)-(2-15) and the modes in [5] 
are identical. The derivation of the modes in [5, section 3.2] is 
a little ambiguous and may be interpreted as only valid in the 
2D case. In reality, the modes (2-14) and (2-15), and thus the 
modes in [5, section 3.2], are the general evanescent plane wave 
solutions of the Maxwell equations in free space, because the 
direction Y defined in Fig. 2 is any direction perpendicular to 
the propagation of the phase. The evanescent plane waves are 
2D in nature (the field is invariant in one direction 
perpendicular to the phase propagation), in the same sense that 
the pure traveling waves are 1D in nature (the field is invariant 
in the two directions perpendicular to the phase propagation). 


II. THE POYNTING VECTOR AND THE VECTORIAL 
RELATIONSHIPS CONNECTING THE £, H, AND k VECTORS 


The Poynting vector gives the flow density of power through 
a surface. It comes from the Maxwell equations [1, 3, 4] and is 
thus valid for any kind of plane waves. Using the physical 
components of the TM mode (2-16) into P=ExH, and 
assuming without loss of generality that XY = 0, the Poynting 
vector is obtained as 
2 42 
Zolto Ay [cosh y (1 + cos 2af)iiy + sinh y sin 2æt ü, | 
(3-1) 
where uy and uy are the unit vectors in directions X and Y. The 
Poynting vector has a transverse component in the direction of 
evanescence Y. This component is alternatively positive and 
negative, and its average value over one period is zero. The 
average value of the X component is % Zo Ho” Ay coshy, which 
is that of a pure traveling wave of magnitude Ho Ay multiplied 
with coshy. The average flow of energy of evanescent plane 
waves is thus in the direction of propagation of the phase_X, as 
with traveling waves, but the instantaneous flow is lying in the 
(X, Y) plane. The TE mode yields the same P vector as (3-1), 


P= 
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with just Eo?/ Zo in place of Zp Ho’. It can be seen that the 
averages on time of (3-1) and its TE counterpart are identical to 
the average Poynting vectors given in [5]. 

Let us now consider the case where the two modes propagate 
and are evanescent in the same directions X and Y. The modes 
are not necessary synchronous which can be expressed by 
replacing their real magnitudes Hp and Ep in (2-14) and (2-15) 
with Ho exp(j ọrm) and Eo exp(j @7z) where ọrm and Orz are 
phases with respect with an arbitrary reference. The physical 
components remain (2-16) and (2-17), with just the argument 
(œt — ox) of the sines and cosines replaced with (œt + ọrm— ox) 
in (2-16) and (œt + ọrg - x) in (2-17). Denoting as (Erm, Hrm) 
and (Ere, Hre) the modified physical fields of the two modes (2- 
16) and (2-17), the Poynting vector 


P= Gy + Ene) x (Hong + Hx) 


can be obtained as: 


(3-2) 


as A? 
P= Z cosh ziz, H; (1 + cos 2(@t + Pry )) 


ES — Ae. 
+ (1+ cos 2(at + Pry)) ux ae 
O 


Es =- 
|Z Hy sin 2(@t + Pry) ) + > sin 2( cat + on) |i : 


o 
+ Ay E, Hy sinh ycosh ysin(gy,—Py,)u, (3-3) 
It can be seen that the X and Y components in (3-3) are the 
addition of the P components of the TM mode (3-1) with those 
of the TE mode, with their respective phases. However, there is 
an additional component in Z direction which vanishes when 
Oru = Ore, 1.e. when the two modes are synchronous. The 
average flow density of energy over one period reads 
= Ae eee ee 
(P) = cosh x) ZH + Bi his 
2 Zi 


T Ay E, H, sinh y cosh ysin(@py, — Ory )U, (3-4) 


Thus, when a wave is composed with the two modes whose 
phases differ with Ag = Oru - Ozz, the average flow of power is 
not parallel to the propagation of the phase, it is lying in the 
plane (X, Z) perpendicular to the direction of evanescence Y. 
Notice that (3-4) can also be derived from the real value of the 
complex Poynting vector 


2 = Real Ën, + Em) x (Hy +Ë) (3-5) 
using the complex fields (2-14) and (2-15) with magnitudes Ho 
exp(/ Ozu) and Eo exp(j Ọre). 

Let us now consider relationships that connect the £, H, and 
k vectors of the two modes. Using wavenumbers (2-11) and the 
TM or TE components (2-14) or (2-15), it can be verified that 
most relationships that hold with traveling waves also hold with 


the TM and TE evanescent modes. For instance, for the TM 
mode, (2-11) and (2-14) yields 


ExB-(2z, cosh? y H, -Ê Z, sinh? z H, |i (3-6) 
c c 


where u, is the unit vector of Z axis. This gives 


k x É =0 m H (3-7) 
Similarly, (2-11) and (2-14) yield 
kxH =-o SE (3-8) 


The same relations also hold with the TE mode (2-15). And 
using (2-11), (2-14), (2-15), it comes 


k-E=0; k.H=0 (3-9) 
for the two modes. The £, H, k vectors of the evanescent modes 
are perpendicular to each other, as the vectors of traveling 
waves in a vacuum. Since £ for the TM mode and H for the TE 
mode are not perpendicular to the propagation of the phase, the 
k vector is not parallel to the propagation of the phase, which is 
also evident from (2-11). The relation of dispersion is obviously 
the same as with traveling waves, i.e. œ = c k, since (2-11) has 
been obtained as a solution of (2-5).This can be verified using 
(2-11): 

2 


k.k=k? = —y (cosh? z— sinh” x) 


Notice that relationships (3-7)-(3-9) can also be derived as in 
[5] from the Maxwell equations, using just the assumption that 
the waveform is (2-1). Other relations identical to those of 
traveling waves can be deduced from (2-11), (2-14), (2-15), as 
the following two for the TM mode and TE mode, respectively: 


Exi=--- Fk 6-10) 
Hoe 


Finally, relationships connecting E£, H and the real part of the 


b u 1 = 
Ex H =——H)k:; 
EO 


wave vector k y (2-11) that governs the propagation of the 
phase can be derived 








pee ye e (3-11) 
E yO E,@ cosh y 
for the TM mode and 
Pes Exk, oe sinh y E,k, (3-12) 
My@ L æ@ cosh y 


for the TE mode. They reduce to (3-7) and (3-8) when y= 0 


(pure traveling waves wherek =k, ). 


IV. EVANESCENT WAVES WITH THE FDTD METHOD 


This section reports numerical experiments that simulate the 
propagation of evanescent plane waves (2-10) by means of the 
full wave FDTD method [7]. To introduce an incident wave in 
a FDTD grid, a useful and widely used method 1s the total- 
field/scattered-field method, also called the Huygens surface 
method [7]. It consists of enforcing electric current Js and 
magnetic current Ks sheets upon a closed surface. Js and Ks 
equal in magnitude the H and F tangential fields that would be 
produced at the surface location if the incident wave were 
present. This permits the incident wave to be reproduced in the 
interior of the surface, whilst no field is produced outside it. 

What is done with pure traveling waves can also be done with 
evanescent plane waves (2-10), for the TM or TE modes (2-14)- 
(2-15). The implementation is similar to that of traveling waves, 
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except the presence of an additional component in the direction 
of propagation, either Ey (2-14a) or Hy(2-15b). Fig. 4-1 depicts 
a 2D domain with a Huygens surface generating the TM mode 
propagating in x direction and evanescent in y direction. In that 
simple case, on the sides C1-C2 and C3-C4 sheet Ks equals the 
transverse component E,, because there is no contribution of Ex 
which is normal to the surface. Conversely, on sides C2-C3 and 
C4-C1, only Ex contributes to Ks. At oblique incidence the 
situation would be a little more complex, the two components 
would contribute on the whole surface. 


The Huygens surface 







e pD 
Y C2 C3 
Incident ieee ee ae 
iT 
TM Mode " eC u 
iT 
i it 
iT iT 
X " iT 
iT iT 
iT iT 
y a} yo nnnnnnnn i ii REE Tí nunnan 20 cells 
i" i 
iT iT 
iT I 
iT "i 
iT iT 
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iT 1" 
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C] Ca 
e A 
X << 
10 cells 


Fig. 4.1. FDTD computational domain for the numerical experiments with 
evanescent waves (Ax = Ay =5 cm, At=0.1 ns). A, B, C, D are the nodes where 
FE, or E, are plotted in Figs. 4.2 and 4.3. The incident wave is evanescent in y 
direction and its phase propagates in x direction. 


For a single frequency incident wave, the components Ex, Ey, 
H, contributing to Ks and Js at (x, y) on the Huygens surface in 
Fig. 4.1 can be expressed as, from (2-16) 


E(x, y,t)=—Z, H, sinh ye? "7° sin(at —x @cosh y/c) 
E (x, y,t) =Z, H, cosh ye” Pees cos(at — x æ cosh zle) 
Hoa e TOR cos(wt — x acosh y/c) 


where yo is the ordinate of the middle of the surface (Fig. 4-1) 
and Ho is the magnitude of H, at y= yo (Ho = 100/Zo A/m in the 
calculations). In the experiment reported in Fig. 4-2, the period 
of the wave was 50 time steps At and its magnitude was growing 
from 0 to Ho over 10 periods to obtain a frequency content as 
close as possible to a single frequency wave. 

The calculated E£, and E, components are plotted in Fig. 4-2, 
at several locations inside and outside the Huygens surface, 
when the steady state has been obtained. The calculation was 
performed with coshy = 1.03, i.e. sinhy = 0.24678. We can note 
that the m/2 advance of E, with respect to E, is well reproduced. 
The magnitude of E, and E, are smaller at point C than at point 
B, because the wave is evanescent in +y direction. The ratios of 
the fields at B and C agree well with the exp[-q/c sinhy (y-yo)] 
evanescence. For example, the FDTD ratio of £, at B and C is 
0.4150 while its theoretical value for the 17 cells separating the 
two nodes is 0.4155. Similarly, the ratio of the FDTD 


magnitudes of E, and E, at point B is 0.2286, in acceptable 
accordance with the theoretical value sinhy/coshy = 0.2396 
(due to the staggered FDTD grid, Ey in Fig. 4-2 1s the average 
of values at the two nodes located Ay/2 from the £, node). 
Outside the Huygens surface, the field leakage is small, as 
illustrated with points A and D. It 1s larger than with pure 
traveling waves and grows with coshy. In this experiment it was 
about 2 % of the field inside the surface. 


Ey Point B 
Ey Point C 
O—0—0 Ex Point B 
o—o—o Ex Point C 
o—o—a Ey Point A 
x—*—x Ey Point D 


E field (V/m) 





80 81 82 83 8 85 8 87 88 89 9 
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Fig. 4-2. Electric field inside and outside a Huygens surface generating a 


sinusoidal TM evanescent field. The locations of points A, B, C, D are shown 
in Fig. 4-1. 


o o o FDTD Results 
Theory 


| Ey (B) 
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Fig. 4-3. £, and £, fields inside the Huygens surface at B and C locations in Fig. 
4.1. The incident waveform at y = yo 1s a Gaussian function. The FDTD results 
(circles) are compared with the theoretical field (continuous lines). 


Beside the above pure sinusoidal wave, it is also possible to 
implement time domain incident waves. Denoting as S(@) the 
spectrum of the desired waveform at a location (xo, yo), the 
spectrum of the field at ordinate y reads, in the case of Fig. 4-1 
where the evanescence is in y direction 


E. (iad O) = 12 H, sinh y S(@) e sinh y (y-yo)/c (4-1a) 
E(x) ,y,@) =Z, H, cosh 7 S(@) e 40%" (4-1b) 
H, (Xp „y, @) =H, S(a) ee sinh yz (y—yo)/c (4-1c) 
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The time domain field at (xo, y) can be obtained by the Inverse 
Fourier transform of (4-1). This field can then be propagated in 
x direction with the phase celerity c/coshy to obtain the 
contributions to Js and Ks at any (x, y) on the Huygens surface. 
Such a time domain incident wave has been implemented with 
a Gaussian pulse in the 2D case 1n Fig. 4-1. 

Results of an experiment with an incident wave of coshy = 
1.03 are reported in Fig. 4-3. The incident pulse at y = yo was 
the Gaussian pulse of spectrum 


l 
S(@) = — =e 
(o) P 

with T= 1 ns. The FDTD fields E, and E, at locations B and C 
in Fig. 4-1 are compared with the theoretical Æ, and E, 
computed by the inverse Fourier transforms of (4-la) and (4- 
1b). As can be observed, the FDTD field inside the Huygens 
surface is in good agreement with the theoretical field. 

Finally, the simple experiments in Figs. 4-1 to 4-3 show that 
evanescent plane waves can be generated in a FDTD grid by 
means of the usual Huygens surface method, and demonstrates 
that such waves can propagate in the FDTD grid in accordance 
with theory. 


-(@T/2) 


(4-2) 


V. EVANESCENT PLANE WAVES IN SOME PROBLEMS 


In many problems of electromagnetism, the analytical 
solution can be interpreted as a set of plane waves. This is well 
documented in textbooks when the solution is propagative, i.e. 
not evanescent. This is the case 1n waveguiding structures 
where the electromagnetic energy of each mode can propagate 
without attenuation above the cutoff frequency. In that regime, 
the field is a set of traveling plane waves. Below the cutoff, the 
modes are evanescent and can also be interpreted as a set of 
plane waves, more precisely a set of evanescent plane waves. 
In the following we address the cases of parallel plate 
waveguides and rectangular waveguides, where the solution is 
a set of two or four evanescent waves (2-14) or (2-15). Another 
example is also addressed, namely the field reflected from 
corrugated surfaces also composed with evanescent plane 
waves. 

The simplest waveguide, depicted in Fig. 5-1, is composed 
with two parallel plates. For the transverse electric (TE) modes 
[3], the £ (x, y, © field is parallel to the plates and in the form 


EzE; sin( "2 5 en (5-1) 
a 


where the time dependence exp(j w?) has been omitted, m is an 
integer, and 


2 2a 
w M 7 


ky =i) -—3 Oe) 
C a 


Below the cutoff angular frequency @.= m n c/a, the quantity 
under the square root is negative and (5-1) becomes: 





22 2 


ma oy 
. (MT A ae 
B,= Eysin{ ™ x) = SS 
a 





(5-3) 


This wave is evanescent in y direction of the waveguide. 
Replacing the sine with exponentials, (5-3) becomes 





(5-4) 


E mn mr a mn mr? a 
J—_~x =y i D She NQ 2 a2 

0 

E, = —-ļe * e!t%*% © ~e 4% e'l S 


Z 2j 


where quantities within the bracket can be viewed as two waves 
Ex and Ez with phases propagating in -x and +x directions, and 
magnitudes decreasing in y direction. Thus, the field E, is the 
addition of two plane waves whose magnitudes are not uniform 
in the plane perpendicular to their propagation. 


= PEC plates a 


Ne 


y 
Fig. 5-1. A parallel plate waveguide. 


Defining now an hyperbolic cosine coshy as 


MTC (5-5) 





cosh y = 
wa 


and using cosh*y— sinh?y=1, the second plane waves in (5-4), 
can be rewritten as 

JE, -j2 coshz x -2 sinhz y 
S e 

2 

This is the waveform (2-10) of a plane wave propagating in 
direction x and evanescent in direction y. Using (5-6) into the 
Maxwell equations and taking account of the invariance in z 
direction, two H components are obtained as 


H „o =-—j/Z, simby E, (5-7a) 
H,, =-1/Z, cosh y E, (5-7b) 


which are parallel and perpendicular, respectively, to the 
direction of propagation x, with same ratios Hy2 / Ez. and Ay, / 
Ex as with the TE mode (2-15). In brief, the second wave in (5- 
4) is the TE mode (2-15). Similar expressions hold for the first 
wave, with just a positive sign in the exponential of the phase 
in (5-6). Finally, below the cutoff frequency, the TE mode in 
the waveguide is a set of two TE evanescent plane waves (2- 
15). 

Let us now address the case of a z-directed rectangular 
waveguide (Fig. 5-2) of section (a, 5) in the plane (x, y). And 
let us consider the component £, of the TE mode (m, n). It is 
known it can be written as 


E (5-6) 





E, =E p sin 27 x] co ZZ ver í (5-8a) 
a b 
where: 
2 22 2D. 
jag a (5-8b) 
Z 2 2 2 
C a b 


Below the cutoff œ. that renders the quantity under the square 
root negative in (5-8b), (5-8a) can be rewritten as 


a| Zt a| 2x7 ) 
m a 7 f a p? pY? 


m a sees 
* Aja 
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ja -—— oy ja gy 
—e í b ) Oz a b E (5-9) 
where: 
ae mn nt o (5-10) 
a’ b? c? 


The four terms in the bracket of (5-9) can be interpreted as four 
plane waves evanescent in z direction, with their phase 
propagating in the (x, y) plane. Consider for instance the fourth 
wave. Its phase propagates in direction O with respect with the 
x axis with 
tan 0 = ZZ (5-11) 
mb 
Denoting as_X this direction of propagation, the fourth wave can 
be rewritten as 
m n 
E, = J E mn m Ne r” eo 8 
4 a 
Similarly, the other three waves in (5-9) propagate in directions 
m0, m+0, and -0. Defining now an hyperbolic cosine as 


(5-12) 


2 2 
TC m n 





cosh y =— | — + — (5-13) 
ie o (a b’ 
the wave (5-12) can be rewritten as 
E -j2coshy X -ĉsinhyz 
A (5-14) 


a 
which is nothing but the waveform (2-10) of an evanescent 
plane wave. 





Fig. 5-2. A rectangular waveguide. 


Concerning the other components of the TE mode, it can be 
seen that £, just differs from £, with factor -n/b in place of m/a 
in (5-9) and (5-14), which means that the resulting £x-y field is 
perpendicular to the direction of propagation @ (5-11). And 
components H, and H, just differ with factors m/a and n/b, 
which means that the resulting H,., 1s parallel to the direction 0. 
Thus, the situation is that of the TE mode (2-15), with a 
magnetic component H,., in the direction of propagation, an 
electric component £x- perpendicular to the propagation and to 
the direction of evanescence z, and with a component H, in the 
direction of evanescence. Concerning the magnitudes of the 
field components, it can be shown that the ratios H-/Fy., and Ay. 
y/Ex-y agree with the ratio of components of (2-15). From which 
the fourth wave in (5-9) is a TE plane wave (2-15). This also 
holds for the other three waves in (5-9). In conclusion, the 
evanescent fields in the waveguide can be expressed as a set of 
four evanescent TE plane waves in the form (2-15) where Y is 


the longitudinal direction z of the waveguide, and X and Z are 
directions lying in the transverse plane. The evanescence 
parameter y (5-13) depends on the sizes of the waveguide (a, b) 
and on the mode (m, n). The same interpretation in term of 
evanescent plane waves also holds with the TM modes. 

Let us now briefly consider another problem where 
evanescent waves are present, the reflection by periodic 
corrugated surfaces addressed in [3]. The reflected field is 
expressed as an infinite sum of Floquet modes. At low 
frequency, when the wavelength is larger than the period of the 
corrugations, the modes are evanescent in the direction x 
perpendicular to the surface. In that case, mode n can be 
expressed in the form [3] 

a in ne z = 2 sin P ‘oy 

H, =H e (2s z) - a 4 g (5-15) 

where 2 b is the period and 0 is the angle of incidence of the 
wave. Defining 





cosh y=sing +~~* (5-16) 
the mode can be rewritten as 
-jÊ coshyz -2 sinhy x 
H,=H,e ° e ° (5-17) 


which is the waveform (2-10). Thus, the Floquet modes are non 
uniform plane waves evanescent in direction x perpendicular to 
the surface and propagating in direction z parallel to the surface. 

In summary, the evanescent fields present in waveguides and 
on corrugated surfaces can be viewed as sets of TM or TE 
evanescent plane wave solutions of Maxwell equations in free 
space (2-14)-(2-15). 


VI. THE SNELL LAW OF THE EVANESCENT WAVES 


In this section and in the next one we address the question of 
the reflection and transmission of an incident evanescent wave 
at an interface between two media. In the above, it has been 
assumed that the medium is a vacuum. However, the results are 
valid as well with a dielectric medium, with just £ in place of £o 
(and the corresponding velocity or impedance in place of c or 
Zo). Let us now consider an interface perpendicular to direction 
z separating two media where c = cı and c = c (Fig. 6-1). The 
incident wave propagates in direction (01, 1), its direction of 
evanescence is nı, and its evanescence coefficient yı. The 
parameters of the reflected wave are denoted as 9,, ,, Nr, Xr, 
and those of the transmitted wave as 02, @2, N2, X2. 

Let us first address the reflection from the interface. The 
phases of the reflected and incident waves must be equal in the 
interface, that 1s ky, = Ay. and ky, = ky. Using (2-12a) and (2- 
12b), we obtain two equations with real and imaginary parts. 
They are equivalent to a set of four real equations for the four 
unknowns 9,, @;, Nr. Xv 


cosh y, cos p, sin 0, = cosh y, cos p, sin 0, (6-la) 
sinh y, cos o, cos 0. siny, +sinh y, sin o, cosy, 

= sinh y, cos ø, cos 6, sin n, +sinh y, sing, cos7, (6-1b) 
cosh y, sin o, sin 0, = cosh 7, sing, sin 0, (6-1c) 
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sinh y, sing, cos@, sin77, —sinh y, cos o, cosy, 
= sinh y, sing, cos 0, sin y, — sinh y, cos p, cosy, (6-1d) 


It can be verified that the following solution satisfies (6-1) 


0,=7-—90, (6-2a) 
P, = Qı (6-2b) 
9, =7-—N, (6-2c) 
t,=-4, (6-2d) 


where (6-2a) and (6-2b) are the same as in the case of a uniform 
wave. The magnitude of the evanescence coefficient y is left 
unchanged, only its sign is changed, which allows the direction 
of evanescence of the incident and reflected waves to be the 
same. 


0 Z 
X2 : Z2 

Transmitted wave Yo 

n2 

y 

X 

Medium 2 Medium 2 
Medium 1 






y 





Reflected wave 


Incident wave Z, x 


Fig. 6-1. The incident, reflected, and transmitted waves at an interface 
perpendicular to direction z. The direction of propagation of the incident wave 
X, is 8; from z axis and its direction of evanescence Yı is nı from y axis (it is 
assumed without loss of generality that ọ = 0). From Fig. 2.3 and (6-2c), the 
direction of evanescence Y,.of the reflected wave is obtained by a rotation 1-1; 
of y axis around_X, axis (y’ in Fig. 2.3 is y here since ọ = 0). 


Let us now consider the transmission through the interface. 
Using (2-12b) and writing that k,2 = 4, in the interface yields a 
complex equation whose real part reads 


cosh y, sing, sin@, =(c,/c,) cosh y, sing, sin8, 


It shows that @2= 0 if ọ1 = 0. Since any situation of the incident 
wave with respect to the interface can be expressed in 
coordinates where ọı = 0, it can be assumed that ©; = @2= 0 
without loss of generality. Writing then that ky = ky and ky = 
k,, and using (2-12), the following three equations are obtained: 


A=cosh y, sin@, =(c,/c,) cosh y, sin 0, (6-3a) 
B =sinh y, cos@, siny, = (c, /c, ) sinh 7, cos 6, sinn, (6-3b) 
C =sinh y, cos7, =(c,/c,) sinh y, cos 7, (6-3c) 


The three unknowns 02, n2, ¥2 can be found from this set of 
three equations. For example, sin@2 and sinhy2 can be found as 





1+D?- (1+ D°} -404C?)4 


sin 0, = z (6-4) 
2(1+C?) 
DANS Ay ACC SBC: 

PES a ee (6-5) 
where D*=A?+B’+C". Then, cosn2 can be obtained from (6-3c) 
as 

cos7, =(c,/c,) sinh y, cos 7, /sinh y, (6-6) 


It can be shown that the quantities under the square roots in (6- 
4) and (6-5) are always positive, so that (6-4) and (6-5) are valid 
with any set c1, c2, 01, N1, ¥1. Equation (6-4) always yields sin@2 
< 1, even in the case of total reflection where 02 = 7/2 (see next 
section). In summary, equations (6-4)-(6-6) are the Snell law of 
the evanescent plane waves. 


VILU. THE FRESNEL COEFFICIENTS 


The reflection and transmission coefficients at an interface 
can be derived as with uniform plane waves, by enforcing the 
continuity of the tangential components of E and H. Here, 
because of the evanescence in direction nı and the presence of 
a component of E or H in the direction of propagation, there are 
four components lying in the interface, even with a pure TM (2- 
14) or TE (2-15) incident mode (except with special values of 
n or 0). This results in four equations of continuity, instead of 
two with uniform TE or TM plane waves. In consequence, the 
continuity cannot be satisfied by a single mode, which means 
that an incident mode yields reflected and transmitted energy 
into the other mode. To derive the reflection and transmission 
coefficients, let us assume the incident wave is the addition of 
the two modes (2-14) and (2-15), with magnitudes Ho; and Eos, 
respectively. And let us denote as Ho, and Eo, the magnitudes of 
the reflected modes, and Ho; and Eo; the magnitudes of the 
transmitted modes. Using then the components of the field (2- 
18)-(2-19) with @1 = @2=0, and taking account of (6-2) which 
in particular means that cos@, = -cos0; and cosn,= - cosn:, the 
continuity of Ex, Ey, Hx, H, in the interface yields the following 
four equations: 


(cos 0, sin 7, cosh y, ~ j sin 0, sinh 7, )Z, (Zs —H,,) 
+ cos @, cos 77, (Ex tele: ) =cos 6, cos77, Ey, 


+ (cos Ø, sin n, cosh y, — j sin 0, sinh y, )Z, H, (71a) 
cos 7, cosh y, Z, Wo; - H,„)-sinn, (E,, + Eo 
=cos7, cosh y, Z, H,, - sinn, Ey, (7-1b) 
(cos 8, siny, cosh y, — j sin 8, sinh y, (Ey - Ea, )/ Z, 

— cos 0, cos 7, (H,, + H,, )=—cos 6, cos7, H,, 
+ (cos 0, siny, cosh y, — jsin 0, sinh y, ) E,,/Z, (7-10) 

sin 7, (H; + H, )+cos 7, cosh y, (E, -E )/ Z, 
=sinn, H,,+cosy,coshy,FE,,/Z,  C-1d) 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


Let us then express the four unknowns Hor, Eor, Hon Eo: of the 
linear system (7-1) in function of the incident magnitudes Ho; 
and Eo; as 


Hy “brt tty Ey (7-2a) 
By Paddy t kaii (7-2a) 
H,, =t,, Ho; +t, Eo; (7-2a) 
Ept Hytta Ej (7-2a) 


The coefficients r and ¢ in (7-2) can then be found by solving 
two times (7-1), one time with Ho;= 1 and Eo;= 0 yields ran, Fhe, 
thn, the, and one time with Ho; = 0 and Eo;= 1 yields ren, ree, ten, 
tee. The following reflection and transmission coefficients are 
obtained 


Tah UZ, Fon -sin 1, 
Fije _ z] W, Z, and l _ -1 V (7-3) 
lah sin 7, lh U,/Z, 
Cie B V l oe W, /Z 
where matrix M and quantities Um, Vm, Wm read 
U,Z, smni U,Z, —sinn, 
M W, Z, -1 W, Z, V, (7-4) 
o |-sinņ, U,/Z, sinn, U,/Z, 
V W, / Z, 7 V, W, /Z, 
U,,= cosh y,,cos7,, 3 V,, = cos@, cos7,, (7-5) 
W = cosh y, cos 8, sin n, — j sinh y, sin 8, (7-6) 


The two solutions (7-3) could be explicitly expressed in 
algebraic form. However, it is simpler to compute them 
numerically, after replacement of the trigonometric and 
hyperbolic functions of 02, %2, n2 with their numerical values 
from the Snell law (6-4)-(6-6). The coefficients plotted in Fig. 
7-2 have been obtained by this method. 

In the following, three special cases where the system is 
simpler are addressed and discussed. 


The special case where the evanescence is in the plane of 
incidence (nı = +90°) 


Let us first consider the special case where nı = 77/2, which 
means that the direction of evanescence Y; is lying in the plane 
of incidence (the plane perpendicular to the interface where the 
direction of propagation X is lying). Equation (6-3c) yields C 
=0 and n2= 7/2. System (7-1) then reduces to 


C, (x1:0,)Z, (Ha - Ho.) =C, (%2:0,)Z2 Ho (1-70) 
Ba T Ly Dy (7-70) 
C(x > 0, ) (Ex; a Eg; \/Z, =C, (xy, > 0- Eo, IL (7-7c) 
add = Fs (7-7d) 

where C\, (7,8) is the special case n = 1/2 of (7-6): 
Cy (x, 0)= cosh y cos 0— j sinh y sin 0 (7-8) 
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System (7-7) is composed with two independent subsets, (7- 
7a)-(7-7d), and (7-7b)-(7-7c), respectively. It means that the 
two modes are decoupled, 1.€. rhe = ren = the = ten = 0. This was 
expected because with the two modes the projections of the £ 
and H fields in the interface reduce to only one E component 
and one H component, with the projections of the TE mode 
perpendicular to those ofthe TM mode. The two modes are thus 
independent. System (7-7) yields the coefficients 


2 C0) 2 Ce) (7-9a) 
Z Cy, 0)+Z, C,(%,4,) 
a LCN) Zi Cas) (7-9b) 


e Z, C,(7,,4,)+Z, C,(7,,6,) 


with in addition tan =1 + rpn and tee =1+ ree. Quantity Cu(X2, 02) 
can be found using (6-4) and (6-5). It can also be obtained by 
means of 
2 1/2 
C AIOS a) (cosh 7v,sin @, + jsinh y, cos 0, y (7-10) 
Ci 
which can be derived using kxı=kx2, with kx from (2-12a). 

In the case 7 =- 7/2, (7-7)-(7-10) remain valid with just the 
change of the sign in front of 7 in (7-8) and (7-10). The 
reflections (7-9) of the cases 7 = 7/2 and 77; =-1/2 are complex 
conjugates. 


b =b =b Q9 aah 


=n 
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Incidence angle 8 (°) 
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Fig. 7-1. Reflection coefficient rj, 1n the case where the direction of 
evanescence is in the plane of incidence (nı = +90°). The interface is between 
a vacuum and a dielectric of .= 4. The case coshy = 1 corresponds to a traveling 
wave of reflection coefficient 7,. There is no conversion of energy from the TM 
mode to the TE mode (rre = 0 for any coshy). 


Reflection ran (7-9a) from the interface between a vacuum 
and a dielectric of €, = 4 (1.e. c1/c2 = Zı/Z2 = 2) is plotted in Fig. 
7-1 for several evanescence coefficients coshy. In the case 
coshy = | the incident wave is a pure traveling wave with the E 
field parallel to the plane of incidence. Its reflection coefficient 
is the r, coefficient of traveling waves in textbooks. 

It can be noticed that the TM mode in the 3D case n= - 17/2 
is nothing but the 2D case considered in a PML medium in [6]. 
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The special case where the evanescence is perpendicular to the 
plane of incidence (n1 = 0) 


Let us now briefly consider the case where n; = 0. Equation 
(6-3b) shows that B=0 and n2 =0. Using ni = 12 = 0 1n (7-1) 
results in a simpler system where the terms with sinn are 
missing, but there remains two equations involving the two 
modes, namely (7-la) and (7-1c). The modes are thus coupled 
at the interface, and the algebraic expressions of the eight 
coefficients (7-2) would remain complicated. As in the general 
case, a numerical solution is simpler. We can just note that in 
this special case the calculation of y2 and Oz is quite simple, 
since (6-3c) yields sinhy2 = C and then from (6-3a) sin@2 = 
A/(1+C’)'? which can also be obtained with B = 0 in (6-4). 

The reflections raa and rre, computed by solving numerically 
system (7-3) are plotted in Fig. 7-2 for a vacuum-dielectric 
interface (e, = 4). For the traveling wave case (coshy = 1) the 
reflection ran is the coefficient r, in textbooks and rre = 0. For 
evanescent waves, part of the energy of the TM incident mode 
is reflected into the TE mode (red coefficients rre in the figure). 


rhh (Coshy = 10) 
rhh (coshy = 3) 
rhh (coshy = 1.5) 
rhh (coshy = 1.1) 
rs= fh (coshy = 1) 
rhe (coshy = 10) 
rhe (coshy = 3) 
rhe (coshy = 1.5) 
o—o—o rhe (coshy = 1.1) 
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Fig. 7-2. Reflection coefficients r, and rre when the direction of evanescence 
is perpendicular to the plane of incidence (nı = 0). The interface is between a 
vacuum and a dielectric of s,= 4. The case coshy = 1 corresponds to a traveling 
wave of reflection coefficient r, (in that case 7,. = 0). For evanescent waves 
(coshy > 1), part of the energy of the TM mode is converted into the TE mode 
(red curves). 


The special case where the reflection is total 


We now consider a case well known in Optics, namely the 
evanescent wave existing behind the interface between two 
dielectric media in the case of total reflection. Let be media 1 
and 2 with speeds cı and c2 such that c2 > cı. And let be a non 
evanescent incident wave, which means yı = 0. The reflected 
wave 1s non evanescent (since y; = -%1 = 0) and equations (6-3) 
become 


cosh y, sin 0, =(c,/c,) sn@, = A (7-11a) 
sinh y, cos 0, siny, =0 (7-11b) 
sinh y, cos y, = 0 (7-11c) 


For 0 such that c2/c; sinO; > 1, (7-11a) can be satisfied only if 
coshy2 > 1, i.e. if y2 # 0. In consequence, (7-11c) implies that 
n2 = + w/2 and then (7-11b) that 02 = + 7/2. Thus, the 
transmitted wave is evanescent with 


cosh vy, =c,/c, sind, (7-12) 


Notice that y2, 02, n2, can also be obtained using (6-4)-(6-6). 
Especially, equation (6-4) gives sin02= 1 since B = C=0 and A 
> | (the inner square root reduces to A?-1). With 02 = n2 = 1/2, 
(2 = 0, and (7-12), the wave vector (2-12) yields ky2 = 0 and 


ko -2 2 sino (7-13a) 
C C 
_ Oo. TET (7-13b) 
ky ==] (c,/c)) sin’ 6, —1 
Cy 
from which the transmitted waveform reads 
-j2 sing x -4 (e2/ 6) sin? 0- z 
w(x,z)=e " e” (7-14) 


This wave propagates in direction x parallel to the interface and 
is evanescent in direction z perpendicular to the interface. The 
reflection coefficient can be found using (7-1). Since the 
incident wave is non evanescent, the angle nı is arbitrary. With 
nı =+ n/2 the degenerated (since sinhyı =0) TM and TE modes 
correspond to E fields parallel and perpendicular, respectively, 
to the plane of incidence. Using nı = 7/2 along with yi; =0 and 
02 = n2 = 7/2 into (7-1) yields two subsets of two equations. The 
two modes are decoupled. The following reflection coefficient 
rın can be deduced as 


Z, cos0, + jZ, 4 (c, /c,) sin? 0, —1 
Z, cos 0, — jZ, 4| (c, /c,} sin? 0, —1 


The same formula applies for ree with just permutation of Z| 
with Z2. As expected, the reflection is total since the modulus 
of (7-15) equals one. Equations (7-14) and (7-15) can be found 
in textbooks on Optics. Using nı =- 7/2, the reflections ran and 
ree Would be the complex conjugates of those obtained with 
m=tn/2. 


(7-15) 


Fyn 


Conservation of energy through the interface 


The average flow density of power that flows from medium 
1 through the interface must equal the average flow density of 
power of the wave transmitted into medium 2. More precisely, 
the < P- > components normal to the interface in media 1 and 2 
must be equal. In general, the transmitted wave is composed 
with the two modes that propagate and are evanescent in the 
same directions. The corresponding flow can be computed 


using (3-4) which yields, in the case of an incident TM mode: 
2 


1 lhe 
(Px) = 5 Cosh 2 H| Z, PA cs 


2 





cos 0, 


+H; taal sinh y, cosh y, sin(A@) cos 7, sin @, (7-16) 








Lhe 


where it is assumed without loss of generality that Ay = 1, and 
where Ag is the difference of the phases of the transmission 
coefficients tp, and the. Concerning the energy flow < P.)> in 
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medium | the situation is more complex since there are three 
modes: the incident wave, either the TM mode or the TE mode, 
and the two reflected modes. Moreover, the incident and 
reflected modes propagate in different directions. Rather than 
deriving an explicit formula that would be complicated, < Pz1 > 
can be easily computed numerically in two steps. Firstly by 
using (2-18) and (2-19) to obtain the x and y complex 
components of the incident (E;, Hi) and reflected (Erm, HrTtm, 
Erre, Hrre) waves. Secondly by using these components to 
obtain < P,; > from the z component of the complex Poynting 
vector as 


(Pa) 


1 SEN = = ae _,* _,* 
S Real| Œ. + Erm + Ere)x(Hi + Hrm +H ) 


Z 


(7-17) 
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Fig. 7-3. Comparison of the average power flow density < P,, > in medium 1 
with the average power flow density < P,. > in medium 2, at an interface 
between two media. The incident wave is the TM mode with evanescence 
coefficient coshy,; = 3. Medium 1 is a vacuum and medium 2 1s a dielectric of 
e, = 4. 


The average power flow densities through the interface 
< P.,> and < P.. >, computed with reflection and transmission 
coefficients from (7-3), are plotted in Fig. 7-3 for an interface 
between a vacuum and a dielectric of €; = 4. The cases nı = 0 
and nı = 90° are the special cases addressed in the above whose 
reflection coefficients are plotted in Figs. 7-2 and 7-1, 
respectively. The case n1 = 45° 1s not special, 1.e. the system (7- 
3) that gives the reflection and transmission coefficients does 
not allow any simplification. As can be seen, in all the cases the 
conservation of energy is well verified, 1.e. the power density 
that flows from medium 1 through the interface equals the 
power density that flows from the interface into medium 2. 


FDTD measure of the reflection coefficient 


The Huygens surface described in section IV has been also 
implemented with an incident wave at oblique incidence 9. In 
that case, the exponential in the spectra of the field components 
(4-1) has to take account of the evanescence in direction 8. This 
experimental setting has been used to compute the reflection 
coefficient of an evanescent wave from an interface between a 
vacuum and a dielectric where the speed of light is c/3. The 
method is the one used with PML interfaces in [8], with a 
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Gaussian pulse incident wave. It consists of using a large 
Huygens surface in y direction, with the interface just behind its 
front side. However, with evanescent waves there appears a 
limitation in the size of the Huygens surface. Consider for 
instance coshy = 1.5, 1.e. sinhy = 1.118, and Ay = 5 cm. For 
frequency 250 MHz, the evanescence coefficient equals 5 10!” 
over 100 cells, and 25 10°° over 500 cells. This means that the 
field is quite inhomogeneous if the FDTD domain is large. 
From this, the numerical noise of the large field region 
propagates toward the region of interest and pollutes it. In 
practice, using a domain a few hundred cells in y direction, it 
has been possible to compute the reflection coefficient up to an 
incidence 8 which depends on coshy, for instance up to 0 = 15° 
with coshy = 1.3. 





Theoretical reflection (7—9a) 
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Fig. 7-4. Comparison of the theoretical reflection of the TM mode (7-9a) with 
the reflection computed with the 2D FDTD method, from an interface between 
a vacuum and a dielectric of speed c/3. For coshy = 1.2 and 1.3 some FDTD 
results which were strongly polluted by the noise due to the non homogeneity 
of the field are not represented (missing circles). 


The 2D FDTD experiments have been performed with the 
TM mode in the case corresponding to the 3D case ni= + 7/2 
whose theoretical reflection is given by (7-9a). Some FDTD 
results are compared with the modulus of (7-9a) in Fig. 7-4 for 
a traveling wave (coshy = 1) and for four evanescent waves 
(coshy = 1.03, 1.1, 1.2, 1.3). An excellent agreement is 
observed, which demonstrates the exactness of (7-9a). 


VIII. CONCLUSION 


This paper is schematically composed with two parts. The 
content of the first part is mainly a revisit of the evanescent 
plane wave solutions of the Maxwell equations in free space. 
They are derived and expressed differently than in the only 
textbook we have found on this question [5]. Especially, 
explicit expressions of the wavenumbers and of the field 
components are provided in Cartesian coordinates which 
renders easy the solution of some problems involving 
evanescent waves. 

The content of the second part is probably more original, 
since we have not found any textbook or paper reporting similar 
content. Firstly the Snell law and the Fresnel coefficients of 
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evanescent plane waves are derived at the interface between 
dielectric media. Secondly, it is shown that evanescent plane 
waves can be enforced in the space domain of such 
computational methods as the FDTD method. This allows 
theoretical results concerning evanescent waves to be verified, 
as in this paper, and opens the way to numerical calculations 
where the source of the electromagnetic energy is composed 
with evanescent waves. 
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